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Abstract 

In this paper, the strong solutions (X, L) of multidimensional stochas- 
tic differential equations with reflecting boundary and possible antic- 
ipating initial random variables is established. The key is to obtain 
some substitution formula for Stratonovich integrals via a uniform 
convergence of the corresponding Riemann sums and to prove conti- 
nuity of functionals of (X,L). 
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1 Introduction and main results 

Let O be a smooth bounded open set in n(x) denotes the cone of unit 
outward normal vectors to dO at x, that is, 

(i) 3 C > 0,Vx g dO^x' g 0,3 k e n(x) 

(x-x',k) + C \x-x'\ 2 >0, (1.1) 
(it) Vx G dO, iSC > 0, 3k G 3f? d , Vx' G 6, 

(x - x', k) + C\x - x'\ 2 > 0, ==>■ k = 6n(x) (1.2) 

for some 9 > , where <90 denotes the boundary of O, O denotes the 
closure of O. We assume that Bt is an valued Tt- Brownian motion 
on a stochastic basis (fl, J 7 , {^"t}te[o,i] ; P) satisfying the usual assumptions. 
We consider the following stochastic differential equations on domain O 
with reflecting boundary conditions: 

X t (x) = x + f b(X s )ds + f a(X s (x)) o dB s - L x t , V t G [0, 1], (1.3) 
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where 6 : 3R d i — 3f? rf and a : i— > x 5R d are continuous functions, o denotes 
the Stratonovich integral. A pair (X t (x),L x ,t G [0,1]) of continuous and 
Tf adapted processes is called a solution to equations (1.3) if there exists 
a measurable set f2 with P(f2) = 1 such that for each uj G Cl 

(i) for each x G (5 the function s i— > with values in 5?^ has bounded 
variation on any interval [0, T] and L x = 0. 

(ii) for all t > 0, X t {x) G 6 and (X t (x),L x ,te [0, 1]) satisfies Eq.(1.3). 
(iii) 

\L x \t= f I {Xa{x )edO)d\L x \ s and L% = [\(X s (x))d\L x \ s (1.4) 
•/ o •/ o 

with £(X s (x)) G n(X s (a;)), where the |T x |t denotes the total variation of 
L x on [0,t]. 

Remark that (iii) implies that the support of d\L x \ is included in {s : 
X s (x) G dO} and the force L x keeps the process X be in O. 

This type of reflected stochastic differential equations has been stud- 
ied notably by Skorohod[17J, Tanaka[20]. Lions and Sznitman[ll], and 
Saisho[l8], and also by Stroock and Varadhan[T9] who used a submartin- 
gale problem formulation, and other authors. Moreover, such reflected 
diffusions can also be reduced to studying multivalued stochastic differen- 
tial equations ( see [H [22l [23] and references therein). It is well-known 
(see [H]) that for any given initial value x G O the Eq.(1.3) has a unique 
solution provided that ||cr(-)|| and |6(-)| are uniformly bounded real-valued 
functions on and satisfy a uniform Lipschitz condition: 3 c > such 
that 

\\a{y) - a(z)\\ < c\y - z\, \b(y) - b(z)\ < c\y - z\ (1.5) 

for any y, z G where bi(x) = 6j(x) + | Yl mt(. x ) a kj(x), \\a 

k,j=i 



Jt{°iM) 2 ^d \b(y)\ := {6,(y)} 2 . 

The natural question aries: does there still exist a pair (Xt,Lt,t G 
[0,1]) of stochastic processes to solve Eq.(1.3) if the initial value is an 
arbitrary random variable Z which belongs to O with probability one and 
may depend on the whole Brownian paths ? 

On one hand, the answer is not immediately clear because one needs to 
deal with anticipating stochastic integration. On the other hand, on a given 
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financial market, different agents generally have different levels of informa- 
tion; besides the public information, some of them may possess privileged 
information, which leads them to make anticipations on some future real- 
izations of functionals of the price process, therefore, for a financial corpo- 
ration, the studying the problem of optimal dynamic risk control/dividends 
distribution has to face the question (see [3 [161 CO E] and references therein). 
The main aim of this paper is to give an affirmative answer to the question 
above. Let us describe now more precisely main results of this paper as 
follows. 

Theorem 1.1. Assume that O is a smooth bounded open set in and 
there exists a function cf> G C\ (5ft d ) such that 

3a > 0, V G dO, VC G n(x), (V0(x),C) < ~aC , (1.6) 

the functions a and b satisfy that a, b and Va are bounded, and the following 

\b(x) -b(y)\ + \\a(x) - a(y)\\ + ||(v<r • a){x) - (Vct • a)(y)\\ 
||(V<7 • V<7 • a){x) - (Va • Vct • o)(y)\\ + ||(v<r -b){x) - (V<7 • b){y)\\ 
+ \\{a T • vV • - (ct t • vV • a)(y)\\ < k\x - y\ (1.7) 

for some constant k > 0, where Co is given by (1.1), a T denotes transpose of 
a, Vcr and V 2 o" denote a's derivatives of first and second order with respect 
to spatial variable x, respectively. Then for any random variable Z with 
P{Z G O} = 1 the pair (X t (Z), Lf ,t G [0, 1]) solves the following stochastic 
differential equation on domain O with reflecting boundary conditions: 

X t {Z) = Z+ fb{X s {Z))ds+ f a(X s (Z))odB s -Lf (1.8) 
Jo Jo 

with X t {Z) G O, and satisfies 

(1) the function s i— ► with values in has bounded variation on any 
interval [0, T] and Lq = 0. 

(2) 

\L Z \t = f I(x si z)edO)d\L Z \ s and L z t = f\{X s {Z))d\L z \ s (1.9) 
Jo Jo 

with £(X S (Z)) G n(X s (Z)), where (X t (x), Lf ,t G [0, 1]) is the unique solu- 
tion of Eq.(l.S), the stochastic integral in Eq.(1.8) is interpreted as antic- 
ipating Stratonovich integral. 

Now we recall the definition of the anticipating Stratonovich integral 
(see |13|). For any t G [0,1], let ir denote an arbitrary partition of the 
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interval [0, t] of the form: tt = {0 = to < ti < ■ ■ ■ < t n = t}. Let 
||7r|| = sup {(tfc+i — ifc)} denote the norm of it. For an ?R d x -Revalued 

0<fc<n-l 

stochastic process / = {f s , s £ [0, 1]}, we define its Riemann sums S-„-(f,t) 
by 

SM,t) = Y,7 T / fsds)-(B tk+1 -B tk ). (1.10) 

We have the following 

Definition 1.1. We say that a stochastic process f = {f s , s € [0,1]} is 
Stratonovich integrable with respect to B if the family S n (f,t) converges in 
probability as \\tt\\ — > 0. In such a case we denote the limit by Jq f s °dB s . 

Let us now describe our approach. To prove Theorem 1.1, the natural 
idea is to replace x in (i), (ii) and (iii) of Eq.(1.3) by the initial random 
variable Z and prove that the pair (X t (Z),L z ) solves the Eq.(1.8). To 
achieve this, the key is to establish the following substitution formula 

a(X s (x)) o dB s \ = f a{X s {Z))odB s , (1.11) 

J o 

f(X s (Z))d\L z \ s = 0, L? = [ £,(X s {Z))d\L z \ s (1.12) 

J o 

for all t G [0,1], where / is a continuous function defined on W 1 with com- 
pact support included in O and £(X S (Z)) G n(X s (Z)). 

The novelty and difficulty of this paper are anticipation, reflection and 
shape of domain O. Since Lions and Sznitman's result in [11] states that 
the solution (X t (x), L%) is Holder continuous of order being less than | 
with respect to the initial value x, the regularity is not good enough to 
satisfy the required hypothesis of substitution formula in the literature 
(see [T3]). it seems that we can not apply the existing substitution for- 
mula to prove (1.11). Moreover, because reflecting boundary conditions 
and shape of domain O, it is also impossible to prove (1.11) by using 
Ito-Ventzell formula used by cone and Pardoux[14J, Kohatsu-Higa and 
Le6n[9]. Instead, we prove (1.11) by showing the uniform convergence 
(w.r.t.x) of the corresponding Riemann Sums S 7T (a(X.(x),t). The Garsia, 
Rodemich and Rumsey's Lemma and moments estimates for one-point and 
two-point motions will play an important role. To prove (1.12) we need 
only to show that the functionals F(t,x) := f(X s (x))d\L x \ s , Lf and 
G(t,x) := J * £(X s (x))d\L x \ s are continuous in (t,x), doing this depends on 
the shape function <p of domain O in (1.6). 
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Remark 1.1. It seems that this new approach can be used to study per- 
turbed stochastic Skorohod equations with anticipating initial processes be- 
cause the solutions of these SDEs are Holder continuous with order be- 
ing less than \ an d not differentiate w.r.t. initial value x (see |7| \ 2 If for 
adapted case). We shall study it in forthcoming paper. 

This paper is organized as follows. Firstly we study the regularity of 
the solution {Xt (x), L x ). Secondly we devote to showing continuity of func- 
tional F(t,x) := Jj f(X s (x))d\L x \ s , Lf and G(t,x) : = jj £(X s (x))d\L x \ s 
. In Section 4 we study moments estimates for one-point and two-point 
motions. In Section 5 we prove the uniform convergence (w.r.t. x) of the 
Riemann Sums S n (a(X.(x),t). Finally we prove Theorem 1.1 in Section 6. 

Throughout this paper we make the following convention: the letter c or 
c(pi,p2,P3, ',•)•) depending only on pi,p2,P3, •, ■, ■ will denote an unimpor- 
tant positive constant, whose values may change from one line to another 
one. 

2 Regularity of the solution (X t (x),L x ) of Eq.(1.3) 

The main aim of this section is to study regularity of the solution (Xt(x),Lf) 
w.r.t. (t, x) via the shape function <p of domain O in (1.6). 

Proposition 2.1. Assume that the smooth bounded open O, the coefficients 
a and b satisfy the same conditions as in Theorem 1.1. (Xt(x), L x ) is a 
solution of Eq.(1.3). Then there is a constant c such that 

E{ sup \X t {x) - X t (y)\ p } < c\x - y|* (2.1) 

0<t<l 

E{ sup m -L v t \ p } < c\x-y\ p (2.2) 

0<t<l 

for any x,y £ O and p > 1. 

Proof. By Holder inequality, we need only to prove Proposition 2.1 for 

p>4. Let bi{x) = bi{x) + \ jr g^xK^x), thatis, 6(x) = &(x) + ±(v<r • 

k,j=i 

o~){x) for any x £ We write solution (Xt(x),L x ) of Eq.(1.3) in Ito's 
form as follows: for x £ O 

X t (x) = x+ [ b{X s )ds+ I a{X s {x))dB s -L x , (2.3) 
Jo Jo 

\L x \t = / I(x s (x)edO)d\L x \ s , (2.4) 
Jo 

L x t = [ £(X s (x))d\L x \ s with C(X s (x)) £ n(X s (x)) (2.5) 
Jo 
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and (Xt(y), L\) also satisfy the same equations above for y € 6. 

Applying Ito's formula to function (f> G C 2 {^ d ) satisfying (1.6) and 
stochastic process X t (x), we have 

4>{X t {x)) = cj){x)+ [\vcj> T a){X s {x))dB s + l\v<i> 1 %){X a {x))d8 
Jo Jo 



(V^ Z)(X 8 (x))d\L x \ 8 

+1 J* tr{ (v 2 0™ T ) (X s (x))}ds, (2.6) 
where tr(A) denote the trace of A. Similarly, we have same expression for 

4{Xt{y)). 

Define f(x) := \x\ p , x = {x\,X2, ■•• , Xd) T € $l d - Then 

V/(x) = p\x\ p - 2 x, y 2 f{x) = p\x\ p - 2 I dxd + p(p - 2)\x\ p ^xx T . (2.7) 

Let m t = X t (x) - X t (y), D t = <j>{X t {x)) + (/>{X t (y)) and N t = exp{-^A}- 
By Ito's formula and (2.5), 

df(m t ) = p\m t r 2 mf(b(X t (x))-b(X t (y)))dt 

+ p\m t r 2 mj (*(X t (x)) - <j(X t {y)))dB t 

- p\m t r 2 m^{X t (x))d\L x \ t 

+ p\m t \ p - 2 mfaXt(y))d\L y \ t 

T{X t {y))){o{X t {x)) - o{X t {y))f 

(2.8) 



+ itr{ V 2 f(m t )(a(X t (x)) - o(X t (y)))(o(X t (x)) - a(X t (y))f}dt, 



dN t = -^N t [(y4> T a)(X t {x)) + (y4> T a)(X t {y))]dB t 
-^N t [(vfb){X t {x)) + (V^iMy^dt 
+^N t (y ( f r i){X t (x))d\L x \ t 

a 

+^N t (v<i ) T o{Xt(y))d\Ly\ t 

a 

-^-N t tr{(v 2 <t>oa T ){X t {x)) + {v 2 <j>aa T ){X t {y))}dt 



+ h Nttr { [F^'KM*)) + (V^a)(X t (y))] T 
x [{vfa){X t (x)) + {vfa){X t (y))] |di (2.9) 
and the stochastic contraction df(mt) ■ dN t is given by 

df(m t ) ■ dN t = -^N t \m t r 2 trl (mf(a(X t (x)) - a{X t (y)))) T 

x [(V<j> T a){X t {x)) + (v4> T v)(X t (y))] T )dt. (2.10) 



Therefore, by Ito's formula again and (2.8)-(2.10), 

N t f(m t ) = exp{-^[0(x)+0(y)]}|x-y| p + f N s df(m s ) 
« Jo 

+ f f(m s )dN s + f df{m s ) ■ dN s 
Jo Jo 

= exp{-^(x)+<f>(y)}}\x-y\P 
a 

+ pf N s \m s \P~ 2 m T s (b{X s {x))-b{X s {y)))ds 
Jo 

+ pf N s \m s r 2 m^(a(X s (x))-a(X s (y)))dB s 
Jo 

- pf N s \m s r 2 (m s ,aXt(x)))d\L*\ s 

Jo 

+ pf N s \m s r 2 (m s ,aMy)))d\Ly\ s 
Jo 

+ I J N s tr{v 2 f(m s )(a(X s (x)) - a(X s (y)))(a(X s (x)) 

- a(X s (y))) T }ds 

- N s f(m s )[(v4> T a)(X s (x)) + (v4> T a)(X s (y))]dB s 

- ^j^N s f(m s )[{V<jH){X s {x)) + (V<^&)(X s (y))]d S 



-x\ 



+ ^j\ a \m 8 \ p - 2 \m 8 \\v<l>(X 8 (x)),^X a (x)))d\L x 

+ |£iV s |m s r 2 |m s | 2 (v0(^(y)),e(^(y)))^l^| s 

- ^ f N s f(m s )tr{{v 2 cf>aa T )(X s (x)) + (v 2 0aa T ) (X s (y))}ds 

2 ft ( 

+ ^2j o N s f(m s )trl[(vc/> T a)(x s (x)) + (vct> T a)(x s (y))] T 

x [(v4> T a)(X s (x)) + (V<f> T a)(X s (y))} jds 

■ j\ s \m s r 2 tr{( m T(a(X s (x)) - a(X s (y)))f 
x [(V</> T a)(X a (x)) + ^ T a){X a (y))] T ^ds 

13 

£ <*(<). ( 2 - n ) 



p2 ,/ 



j=i 
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By condition (1.6), 



^|m s | 2 (v0(X s (x)),£(X s (x))) - (m a ,£(X t (x))) < 0, d\L x \ s a.s. 
^|m s | 2 (v0(X s (y)),e(X s (y))) + (m s ,t(X t (y))) < 0, d\V\ 8 a.s. 



Hence 

o 4 (t) +og(t) < 0, o 5 (*)+oi (0<0. (2.12) 
Using G C 2 (K d ) and (2.11)-(2.12), 

(m*(t))* < cE«(i)) 2 + J>*(i)) 2 + £>*(*)) 2 ], (2.13) 

i=l i=6 i=ll 

where a*(t) = sup se[(y] {|a;(t)|}, m*(t) = sup se[ot] {|m(i)|}. 

Since 4> is bounded, by Burkholder inequality(see [2]) and (1.7), we have 

E{(a^)) 2 } < cEyy*\m s \^-hr{lmT(a(X s (x))-a(X s (y)))f 

x[mJ(a(X s (x))-a(X s (y))))}ds 

< cE j jT Im.l^UKX.^)) - a(X s (y)))|| 2 d S 

< c rE{|m*(s)| 2p }(i S . (2.14) 

JO 

Similarly, since <j> and V0 a are bounded on (5, we also have 

E{(4(£)) 2 } < c rE{|m*(s)| 2p }ds. (2.15) 
J o 

Using <j) and a are bounded on O, the condition (1.7) and Holder in- 
equality, 

E{(a? 3 (i)) 2 } < cE|^iV 2 K| 2 ^ 4 |tr|KKX s (x))-a(X s ( y )))) T 
x[(v^)(X s (x)) + (V0 T a)(X s (y))] T j| 2 d S 

< ce| jf |m,| 2p - 2 ||(a(X s (x)) - a(X s (y)))fds\ 

< c f V{\m*{s)\ 2rp }ds. (2.16) 

J o 



By the same way as in (2.16) 

E{«(t)) 2 } < cfv{\m*{s)\*}ds (2.17) 
J o 

for i = 2,6,8,11,12. 

Putting the above inequalities (2.13)-(2.17)together implies that 

E{(m*(t)) 2p } < c\x-y\ 2p + c B{(m* (s)) 2p }ds, (2.18) 

Jo 

By Gronwall inequality, 

E{ sup \X t (x) - X t (y)\ 2p } < c\x-y\ 2p . (2.19) 
*e[o,i] 

So proof of (2.1) has been done by Holder inequality. Using 

\L% - L\\ <\x-y\ + \X t (x) - X t (y)\ + | f (b{X s {x)) -b(X s (y)))ds\ 

Jo 

+ | / (a(X s (x)) - a(X s (y)))dB s \, 
Jo 

(2.2) is a direct consequence of (2.1). Thus we complete proof of Proposition 
2.1. □ 

Proposition 2.2. Assume that the smooth bounded open O, the coefficients 
a and b satisfy the same conditions as in Theorem 1.1. (X t (x),Lf) is a 
solution of Eq.(l.S). Then there is a constant c, which is independent of x, 
such that for any p > 2 

supE{|^(x)-X s (a;)| 2p } < c\t - s\% , (2.20) 
sup_E{\Lf - L x s \ 2p } < c\t - s\i . (2.21) 

Proof. By Holder inequality, we need only to prove Proposition 2.2 for 
p > 4. For t > s > 0, similar to that of Proposition 2.1, we define mt,Dt, 
Nt and / here by 

m t (x) = X t (x)-X s (x) 

t_ rt rt 

b{X u (x))du+ / a(X u (x))dB u - / aX u (x))d\L x \ u 

J s J s 

with£(X s (x)) G n(X s {x)), 

A = <f>(X t (x)), 

N t = exp{-^A}, 
a 

G t = N~\ 
f(x) = \x\ 2 ,x = (xi,x 2 ,--- ,x d ) T G?R. d . 
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By the same way as in (2.11), 

rt 

\X t {x) - X s {x)\ 2 = 2G t N u m^b{X u {x))du 

J s 

+ 2G t I N u mla{X u {x))dB u 

J s 

- 2G t ! N u (m u ,Z(X u {x)))d\L x \ u 

J s 

+ G t f N u tv{{aa T ){X u {x))}du 

J s 

- -G t f N u f(m u )[(v4> T a)(X u (x))}dB u 

2 

a 



-G t I iY u /(m u )(V(/> T 6)(X u (x))du 



+ -G t I N u \m u \ 2 (V(l>(X u (x)),{;{X u (x)))d\L x \ u 

J s 

- -G t / f iV u /(m u )tr{(v 2 0(Ta T )(X u (x))}dn 

^ J s 

+ ^GtJ s N u f(m u )tv^[(V(f> T a)(X u (x))} T [(V<p T a)(X u (xy 

- -G t f N u tr{ (ml(a(X u (x)))) T ((V^a^Xuix^fldu 



a 
10 



£*(*)■ 



By condition (1.6), 

d 3 (t) +d 7 (t) < 0. (2.22) 

Therefore 



B{\X t (x) - X s {x)\^} < c(p)^{m)\ P } + c(p)J2E{\d i (t)\ P } 

i=l i=4 
10 

+c(p)J2v{\di(t)\ P }- (2-23) 



i=7 



Since a, Nt and Gt are uniformly bounded, by Burkholder (see [2j) and 
Holder inequalities and Young's inequality: for any real positive x, y, rj,p, q 
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with p 1 + q 1 = 1 there exists c < +oo such that xy < r/x p + cy q , we have 
E{|d 2 (t)| p } < cE{| f N u m T u a{X u (x))dB u \ p } 

J s 

< cE{ / ^tr{a T (X M (x))m u m^(X„(x))}d«}^ 

./s 

ft p 

< cE{ / \m u \ 2 du} 2 

J s 

< c\t-s\2+cf E{\m u \ 2p }du. (2.24) 

is 

Similarly, 

E{|d 5 (i)| p } < c f E{\m u \ 2p }du. (2.25) 

is 

Since <r, V0cr, iVf and Gt are uniformly bounded, by Holder inequalities 
and Young's inequality, we have 



E{|d 10 (t)| p } < cEUj ^tr|K(a(I tl (x)))) I ((V^)(I u (i))) J \du\ p 

< ce| ^ |m u |dul 

< c\t -s\ p + cj^ E{\m u \ 2p }du. (2.26) 

Similarly, 

E{|(ii(t)| p } < c|t-s| p + c J E{\m u \ 2p }du, (2.27) 

E{\d±(t)\ p } < c\t-s\ p , (2.28) 

E{\di(t)\ p } < c I E{\m u \ 2p }du, i = 6,8,9. (2.29) 

Putting the above inequalities (2.23)-(2.29) together, we obtain 
E{\X t {x) - X s {x)\ 2p } < c\t - a|§ + c / E{\X u (x) - X u {x)\ 2p }du. 

J s 

The Gronwall-Bellman inequality (see |15j for Theorem 1.3.1) implies that 

E{\X t {x) - X s (x)\ 2p } < c\t - s\2 . 
Therefore the proof of (2.20) has been done. Using 

\I%-L*\ < \X t (x) - X s (x)\ + | fb(X u (x))du\ + | ! a(X u (x))dB u \, 

J s J s 
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(2.21) is a direct consequence of (2.20). Thus we complete proof of Propo- 
sition 2.2. □ 

Since the domain O is bounded, the following follows immediately from 
Proposition 2.1 and Holder's inequality. 

Proposition 2.3. Assume that the smooth bounded open O, the coefficients 
a and b satisfy the same conditions as in Theorem 1.1. (Xt(x), Lf) is a 
solution of Eq.(l.S). Then there is a constant c, which is independent of x, 
such that 

E{ sup \X t (x)\ p } <c(l + \x\) p , (2.30) 

0<t<l 

E{ sup \L X \ P \ < c(l + \x\) p (2.31) 
0<£<1 

for any x € O and p > 1. 

3 Continuity of functionals of local times 

Proposition 3.1. Assume that the smooth bounded open O, the coefficients 
a and b satisfy the same conditions as in Theorem 1.1. (Xt(x), Lf) is a 
solution of Eq.(l.S). Then the functions Xt{x) , Lf, F(t,x) and G(t, x) are 
jointly continuous in (t, x) on [0, l]xO, where F{t,x) := Jj f{X s (x))d\L x \ s 
and G(t,x) := J* * ^(X s (x))d\L x \ s , f is a continuous function defined on 
with compact support included in O and ^{X s (x)) G n(X s (x)). 

Proof. By Kolmogorov's continuity criterion of random fields (see [10] 
for Theorem 1.4.1 ), Proposition 2.1 and 2.2, the functions Xt(x) and Lf are 
Holder continuous in (t,x). Since proof of continuity of G(t,x) w.r.t.(t, x) 
is similar to that of G(t,x), we need only deal with the proof of F(t,x). 
Remarking that 

\F(t,x) - F(s,x)\ < sup{|/(y)|}|L? - L% (3.1) 

the function F(t, x) is continuous in t uniformly with respect to x in com- 
pact set O by Proposition 2.2 and Kolmogorov's continuity criterion( see 
Theorem 1.4.1 in [ID] ). Thus, it suffices to show the continuity of F(t,x) 
w.r.t.x for any fixed t. Let x n ,x £ O with x n — ► x as n — ► +cxd. By 
Propositions 2.1-2.2, and Xt(x) and L x are Holder continuous in (t,x) € 
[0, 1] x 6, we have 
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uniformly in t, as n — > +00. Therefore, there exist constants C\, C > 1 
such that for all n > 1 

\L Xn \ < C+ \L X \ < C+ \L X \! < C + d (3.3) 

due to bound of total variation of L x on [0, 1]. Since the function f(x) is 
bounded and continuous, by (3.2) and (3.3), 

r \f(X s (x n )) - f(X a (x))]dL* a »\ — > (3.4) 



/ 

J 



as n — > +00. Because L Xn and L x are continuous processes with bounded 
variation, by (3.2), the sequence of finite sign measures dL Xn on [0,1] 
converges weakly to the finite sign measure dL x on [0,1]. Therefore, for 
bounded continuous function f(X s (x)) on [0, 1], we have 

hm / f(X s (x))dL x " = [ f(X s (x))dL x . (3.5) 
n ^°°Jo Jo 

The proof of Proposition 3.1 follows from (3.4) and (3.5). □ 

As a direct consequence of Proposition 3.1, we have the following. 

Proposition 3.2. Assume that the smooth bounded open O, the coefficients 
a and b satisfy the same conditions as in Theorem 1.1. (X t (x),Lf) is a 
solution of Eq.(l.S). Then there exists a set fl £ T with P(fi) = 1 such 
that for each uj € $7 

f b{X s (x))ds\ x=z = f b(X s (Z))ds, (3.6) 
Jo Jo 

f f(X s (Z))d\L z \ s = 0, Lf= j\{X s {Z))d\L z \ s (3.7) 
Jo Jo 

for all t G [0,1], where f is a continuous function defined on %l d with 
compact support included in O and £(X S (Z)) £ n(X s (Z)). 

4 Moments estimates for one-point and two-point motions 

For any R > and x <G [-R, R] d D O, let (X t (x),Lf) be a solution of 
Eq.(1.3). We define S v (t,x) and I(t,x) by 

S n (t,x) := S n (a(X.(x)),t), 



I(t,x) := [ a(X s (x))odB s 
Jo 

= [ a{X s {x))dB s + \ I {Va-a)(X s {x))ds. 
Jo 1 Jo 
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Write 



n-l 
fc=0 

+ Y, tk+1 _ tk (J t (v(Xs(x))-a(X tk (x)))ds)(B tk+1 -B tk ). 

(4.1) 

By Ito's formula and (1.4), for s > tfc, 

(Jij(X s (x)) - (jy(X tfc (x)) = / (Very • a){X u (x))dB u 

+ / (\/a ir b){X u {x))du 
~ f ■ *i){X u {x))d\L x \ u 

\J t {(v^-(Va-a))(X tt (a;))}du 

So we informally write a(X s (x)) — a(X tk (x)) as follows: 
a(X s (x)) - a(X tk (x)) = f (Va ■ a) (X u (x))eZ£ u 

+ / (v<r-6)(X u (ar))du 

■>t fc 

+ |(v<7-(V<7-<r))(X u (a;))|d« 
+ ijf trj(vV-aa T )(X u (x))jdu. 

(4.2) 

Thus we can write S^i, x) — x) as follows: 



+ 
+ 
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S„(t,x)-I(t,x)=J2 A ™, ( 4 - 3 ) 



i=i 
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where 
A l7r (x) 



J2<j(X ti (x))(B ti+1 -B ti ) - f a(X s (x))dB s , 



=0 
n-l 



=o 

1 f 

2 Jo 

n-l 



(Vcr • a)(X s (x))ds, 
I rU+i 



=o 

n-l 



I fU+i 



i ^ i f u+i ( r *\ 

2T, t . +1 _ t J t ds \j t {y<T-y<r-a)(X u (x))duj(B ti+1 -B ti ) 



i=0 

n-l 



Atoix) := IJ^J^tJ^ 1 ds{j\v{(v 2 a-a-a T ){X u {x))}d u yB ti+1 -B ti 

Proposition 4.1. Assume that the smooth bounded open O, the coefficients 
a and b satisfy the same conditions as in Theorem 1.1. (X t (x), Lf) is a 
solution of Eq.(l.S). Then for any p > 2 and R > there exist constant 
c(p,R), which is independent of t and n, and (3o € (0, 1) such that 

sup B{\S 7r (t,x)-I(fx)\ 2 P}<c(p,R)\\7rfoP. (4.4) 

xe[-R,R] d no 

Proof. By Burkholder-Davis-Gundy and Holder inequalities, we have 



{E{\A l7T (x)\ 2 P}}? < c(p) 



n-l ,t i+1 \ p-\ 



E / \(o-(X s (x))-a(X k (x))\ 2 ds 



=o ti 



n-l 



i=0 
n-l 



< c(p)J2 El / |(cr(X s (x)) -a(X ti (a;))| 2 dsP 



n— 1 r /<t i+1 
r\ L •/ ti 



i=0 

< c(p)||7r||5. 



where we have used Proposition 2.2 and the condition (1.7). Thus 

E{|^ l7r (z)| 2 n<c||7r||!. (4.5) 
Using Fubini Theorem, A-^ can be further written as 



A 27T (x) = A%(x)+A%(x)+A%(x) 



(3), 



(4.6) 
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where 

™ _1 i rU\i 

I ; i l — / , / . 



450*0 == Er^ZT l*\u+i-u){{Va-a){X u {x))-{Va-a){X u {x)))du, 
4?0*0 := -jE/ ((^.a)(X u (x))-(Va.a)(X ti (x)))du, 

n — 1 s ^ fti+l 

4?(*) := E{^-^— tSl ^- u )(y°-°K X u{x))dB u ){B ti+1 -B ti ) 

1 Z"*^ 1 1 

r( / (U+i - u){Va ■ a)(X u (x))du) \ 

H+l — H hi J 

n-1 

:= ^(Ai + Bi). 



i=0 

It follows from (1.7) and Proposition 2.2 that 



{E{l4JWr}}' < E / l+1 {E{||((v f T.a)(X tt (x))-(Va.a)(X ti (x)))n}i^ 

i=0 * l 
n—1 .(. j 

< c£ / {E{\\X u (x)-X ti (xW}}ku 



< c||vr|| 4 . (4.7) 

Similar arguments lead to 

{E{|4 2 JWI P }}"<c||vr||i. (4.8) 

Since 45 is a martingale and Vo" • a is bounded on (5, using Burkholder- 
Davis-Gundy and Holder inequalities, we obtain that 

n— 1 



{ E {i45(*)i 2p }}"< c E{^iA+^i 2 r 

^ i=0 
n— 1 n—1 

< c^(E|^| 2 ^ +c^(E|5 i | 2 P)p 

i=0 i=0 

= c(p)£ (E{|— — - / (t m -u)(V<r-a)(X u (x))d£ u ) 

+ (E{ | ^-i-^- ( (i m - u) ( Va • a) (X u (z) )du \ }^) F } 

< c(p)||7r||. (4.9) 
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So we deduce from (4.7)-(4.9) that 

E{|A 27r (^)| 2p }< C ||7r||i (4.10) 

By Propositions 2.1 and 2.2, it follows from Kolmogorov's continuity crite- 
rion( see Theorem 1.4.1 in [10] ) that there exist a random variable K with 
E|i^(u;)| p < +oo and a positive constant G (0, 1) such that 

sup \L x t -L x s \<K{u)\t-sf. (4.11) 
xe[-R,R] d no 



B{\A 37T (x)\ 2p } < E{(V|- -/ is / (V(T.0(^«(a;)Wlur 

^ ti+i - U J t . J t . 

71-1 

x(^is t!+1 -^i 2 n 

i=0 

n— 1 Ti—l 

< C E{(^(i^| ti+1 -i^yy x(^|B tj+1 - jBti iy} 



8=0 1=0 

71-1 



< cE{(sup|L£ +1 -Ll\\L*\ x ) p x (^|B ti+1 -B ti | 2 ) P } 

i ■ „ 

4=0 

< c(E{(sup|2£ +1 -^|) 3p })§(E{|^|f})l 



71-1 



(E{(^|^ +1 -^| 2 )n) s 



i=0 



< cjlvrfP (4.12) 

due to the inequality (4.11) and Vcr is bounded on O. For p > 1, by Holder 
inequality and S7 2 a ■ a ■ a T is bounded on (5, 



< 



E{|^(x)| 2 ^} 

^(Elr^ r +1 ^ rTr{(vV.a-^)(X u (x))}^| 2 ) 



i=0 
n-1 

- ' 

/_() 

71-1 



< c||7rfsupE{(^|^ +1 -i? ti | 2 ) p } 

i=0 

< c||vrf. (4.13) 
Similarly, 

E{\A i7T {x)\ 2p } < c\\ir\\ p , fori = 4, 5. (4.14) 
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Putting the above estimates (4.5), (4.10) and (4.12)-(4.14) for A in (i = 
1, • • • ,6) together, we deduce that 

sup E{\S n (t,x)-I(t,x)\ 2 P}<c(p,R)\\7r\f°P, 

xe[-R,R] d no 

where (3q = min{^,/3}. The proof of Proposition 4.1 is complete. □ 

Next result is the moment estimates for the two point motions. 

Proposition 4.2. Assume that the smooth bounded open O, the coefficients 
a and b satisfy the same conditions as in Theorem 1.1. (Xt(x),Lf) is a 
solution of Eq.(l.S). Then for any p > 2 and R > there exists constant 
c(p,R), which is independent oft and n, such that 

E{ sup \S n (t,x)-S n (t,y)\P}<c(p,R)\x-y\P, (4.15) 
te[o,i] 

for all x, ye [-R,R] d DO. 
Proof. Similarly as (4.3), 

6 

S*(t, x) - S w (t, y) = J2 A -( x > V)> ( 4 - 16 ) 

i=\ 

where 

n-1 



A ln (x,y) := H X tM)) ~ <X u {y))] {B u+1 - B u ), 

i=0 

A 2w (x,y) := _ f . / ds \ / [( V 

a • a) (X u (x)) 

i=0 i J t% ^ Jti 

-{Va ■ a)(X u (y))]dB u \(B ti+1 - B ti ), 

"--1 1 pt i+ l ( fS 

A 3 n(x,y) := V- -/ ^ / (Va-O(^(^)^l^k 

- £ (7a • (X u (y))d|Lf| u | (B ti+1 - B u ) , 

1 f*^ 1 f /" s 
Ator(x,y) := V- -/ ds\ / [(Va • b) (X u (x)) 

-(Va • b)(X u (y))]du\{B u+1 - B u ), 
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1 2— i 1 z"*^ 1 f f s 

A 57r (x,y) := dsl [(Va ■ Va • a) (X u (x)) 

-(Va • Va • <7)(X u (i/))]du - Bt.), 

:= ±£_L_ r +1 dJ f S Tr{(v 2 a.a-a T )(X u (x)) 
-(v 2 a • a • a T )(X u (y))}dn|(B ti+1 - B u ). 

Let \L X — L y \ t denote the total variation of L x — L y on [0, t] for any t € [0,1], 
by (1.4), we have for any s € [<i,ti+i] 

| f {Va-t){X u (x))d\L x \ u - f \va-t){X u (y))d\Ly\ u \ 2 

S {Va)(X u (x))dLZ - f (Va)(X u (y))dLl\ 2 

i J to 



< c\J [(Va)(X u (x))-(Va)(X u (y))]dL* u \ z + \j (Va)(X u (y))d(Ll - Ll)\ z 

< c sup {\X s (x) - X s (y)\}{\L x \t i+1 ~ \L*\ ti } + c(\L* - #Vi - \L* - V\ u }. 

se[o,i] 



So 



n-l 



i=0 i+1 ~~ ' Jt * Jt i 

. s n-l 

- (V«7 ■ e)(X tt (y))d|L"|2}tfa x ^ |B ti+1 - B t 

ti „- n 



i=0 



< 



2c 2 snp {\X s {x) - X s {y)\ 2 }\L*\\\ 
se[o,i] 

+2c 2 sup {iLf-^IKlLHi + l^li) 



S6[0,l] 



n-l 



i=0 
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It follows from Propositions 2.1-2.3 and Holder's inequality that 

E{ sup \A 37T (x,y)\ 2 P} 
te[o,i] 



< c(E{ sup \X a (x) - X s (y)f P })H^{\L x \t P }^ 

0<s<l 
n-l 

x(E{(El^ +1 -^| 2 ) 3P }) 5 



i=0 



+c(E{( sup |LJ-XJ|*})t(E{|L*|*+|L«'|*})i 



0<s<l 
n-l 



:(e{(Ei^i-^iT}) 



i=0 

< c|a;-y| p . (4.17) 

By Burkholder-Davis-Gundy inequalities, the condition (1.7) and Proposi- 
tion 2.1, it follows easily that 

E{ sup |A l7r (x,y)r} <C(p,R)\x-y\ p . (4.18) 
te[o,i] 

Using Burkholder-Davis-Gundy and Holder's inequalities, the condition 
(1.7) and Proposition 2.1, 



\ — Til /*t f P s 

E{ sup \A 27T (x,y)\P} Y / ds\E{\ / [(Va ■ a) (X u (x)) 

te[o,i] J i=0 h+i ~ ^ Jti I 



-(Va-<7)(X tt (y))]dB Tt j(B ti+1 - J B ti )| 1 ' 

< Eir^TT f 1+1 dsU{\ f S [(Va.a)(X u (x)) 

J_ J_ 

-(Va • *)(X u (y))]dB u \ 2 ^ 2P (b{\(B u+1 - B u )\^ 2P } 

< Eir^-r f tl+1(ls ( f {n\Xu( X )) - x u (y)\^}}hu\ 

x(t i+ i 

< c(p)|x-y|. (4.19) 
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Similarly, 

(E{ sup \A 4n (x,y)+A 57r (x,y)\P}) P < c{p)\x - y\, (4.20) 

V te[o,i] / 

i 

(b{ sup \A e7r (x,y)\v}Y <c(p)\x-y\. (4.21) 

V te[o,i] / 

Thus, combining above estimates (4.16)-(4.21) together , we complete the 
proof. □ 

The following result can be proved similarly as Proposition 4.2, but its 
proof is very easy, we omit it here. 

Proposition 4.3. Assume that the smooth bounded open O, the coefficients 
a and b satisfy the same conditions as in Theorem 1.1. (Xt(x), Lf) is a 
solution of Eq.(l.S). Then for any p > 2 and R > there exists constant 
c(p,R), which is independent oft and n, such that 

E{ sup \I(t,x)-I(t,y)\P} <C(p,R)\x-y\P, (4.22) 
te[o,i] 

for allx,y<= [-R,R] d nO. 



5 Uniform convergence of the Riemann sums 



Let R > and p > 1 be given. Define G R := [-R, R] d n O. Then the 
following is a direct consequence of Garsia-Rodemich and Rumsey's Lemma 
(ef.pjUEllSJ). 

Lemma 5.1. Let f : £1 x ffl 11 — > K n be a measurable stochastic field taking 
values in ffl 1 which is continuous, P- a.s., p > 1. Then there exists a 
constant c(p, R) such that 



E{ sup p{f{x)J(y)Y) 

x,y£G R 



<c(p,R) 



E f p(f(x),f(y)\y 
d(x,y)P 



I {x -i y} dxdy, (5.1) 



>G R xG R 

where (3? m ,ci) and (5R n ,p) are metric spaces. 

Now we prove the main result of this section. 

Theorem 5.1. Assume that the smooth bounded open O, the coefficients 
a and b satisfy the same conditions as in Theorem 1.1. (X t (x), L^) is a 
solution of Eq.(1.3). Then for any p>2 and R > 0, 

lim E{ sup \SJt,x) -I{t,x)\ 2p \ 



0. 



(5.2) 
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Proof. Since sup xeGa {\f{x)\} < aup XtyeGH {\f{x) - f(y)\} + |/(x )| for 
any xo € Gr and function / on we have 

E{ sup \S v (t,x) -I(t,x)\ 2p } 

<c(p)E{ sup \S n (t,x)-S n (t,y)-I(t,x)+I(t,y)\ 2 P} 

x,yeG R 

+c(p) sup E{|^(t,x)-/(t,x)| 2 f} 

= B l7r + B 2w . (5.3) 
By Lemma 5.1, 

(5.4) 

By Propositions 4.1, 

f |^(t, x) - S n (t, y) - I(t, x) + I(t, ) c^rf°P 

E \ j 7 ^> - -> ° 

as ||7r|| — > 0. 

Therefore, by Propositions 4.2, dominated convergence theorem and (5.4), 
we have 

B l7T — ► , as 7T -> 0. (5.5) 

By Proposition 4.1, 

B 2 tt — ► , as 7T -> 0. (5.6) 

Thus we complete the proof by (5. 3), (5. 5) and (5.6). □ 

6 Proof of Theorem 1.1 

We will prove that (X t (Z),Lf ) solves the anticipating reflected SDE (1.8). 
Since (X t (x),Lf) is a solution of Eq.(1.3), by Proposition 3.1, we need 
only to prove (1.11). By Theorem 5.1, the following holds almost surely on 
{uj;Z(uj) <E G M } 

ft 



j a(X s (x)) o dB a \ zX { U ;Z{ u )eG M y 
Jo 



limS n (t,x)\ x=zX{ ^z^)eG M } 
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= f a(X s {Z))odB a X{ u . Z {u>)ea M }- 
Jo 

Letting M — > oo, we obtain the substitution formula (1.11), and therefore 
prove the Theorem. □ 

Remark 6.1. By checking carefully the proof of Theorem 1.1 and using 
Theorem 3.1 proved by Lions and Sznitman (see II If ), the conditions on 
O in Theorem 1.1 can be weaken, that is, if O satisfies the admissibility 
condition (see \llf . page 521) and the following condition: there exists a 
function <j) in Cf(5R d ) such that 3a > 0, Vx G dO, Vy G O, V£ € n(sc) 
=^ a(y < H a:; )> 0\v ~ x \ 2 ~ (v ~ x >0 — Theorem 1.1 also holds. 



Acknowledgements. This work is supported by NSFC and SRF for 
ROCS, SEM. The author would like to thank both for their generous fi- 
nancial support. 

References 

[1] Asmussen, A., H^jgaard and Taksar. M.: Optimal risk control and 
dividend distribution polices, example of excess-of loss reinsurance for 
an insurance corporation. Finance and Stochastics, 4(2000)299-324. 

[2] Barlow, M. T.; Yor, M.: Semimartingale inequalities via the Garsia- 
Rodemich-Rumsey lemma, and applications to local times. J. Funct. 
Anal., 49, 198-229(1982). 

[3] Baudoin, F. and Nguyen-Ngoc, L.: The Financial value of weak infor- 
mation on a financial market. Finance and Stochastics, 8(2004)415-435. 

[4] Cepa,E. and Lepingle, D.: Diffusion particles with electrostatic repul- 
sion. Probability Theory and Related Fields, 107 (1997)429-449. 

[5] Cepa,E. : Probleme de skorohod multivoque. Ann. Probab. 
26(1998)500-532. 



23 



[6] Cadenillas, A., Choulli, T., Taksar, M. and Zhang, L. : Classical and 
impulse stochastic control for the optimization of the dividend and risk 
polices of an insurance firm. Mathematical Finance, 16(2006) 181-202. 

[7] Doney, R.A. and Zhang, T.: Perturbed Skorohod equations 
and perturbed reflected diffusion processes. Ann.I.H. Poincare-PR, 
41(2005)107-121. 

[8] Imkeller, P. and Scheutzow, M. : On the spatial asymptotic behavior 
of stochastic flows in euclidean space. Annals of Probability,27, 109- 
129(1999). 

[9] Kohatsu-Higa, A.; Leon, J. A.: Anticipating Stochastic Differential 
Equations of stratonovich Type, 36(1997)263-289. 

[10] Kunita, H.: Stochastic flows and stochastic differential equations. 
Cambridge University Press, 1990. 

[11] Lions, P.L. and Sznitman, A.S. : Stochastic differential equations with 
reflecting boundary conditions. Comm. Pure Appl. Math. 37(1984), 
511-537. 

[12] Liang, Zongxia : Spatial asmptotic behavior of homeomor- 
phic global flows cor non-Lipschitz SDEs. Bull. Sci. math. (2007), 
doi: 10. 1016/j.bulsci.2006. 12.001. (In Press). 

[13] Nualart, David : The Malliavin calculus and related topics. Probability 
and its Applications (New York). Springer- Verlag, New York, 1995. 
xii+266 pp. ISBN: 0-387-94432-X. 

[14] Ocone, D and Pardoux,E: A generalized Ito-Ventzell formula. Appli- 
cation to a class of anticipating stochastic differential equations . Ann. 
Inst. H. Poincare Probab. Statist. 25:1 (1989) 39-71. 



24 



[15] Pachpatte, B.G.: Inequalities for Differential and Integral Equations. 
Academic Press, 1998.ISBN0-12-543430-8 

[16] Pikovsky, I. and Karatzas, I.: Anticipative portfolio optimization. Ad- 
vances in Applied Probability, 28(1996)1095-1122. 

[17] Skorohod, A.V.: Stochastic equations for diffusion processes in a 
bounded region I. Theory Probab. Appl. 6 (1961)264-274. 

[18] Saisho, Y.: Stochastic differential equations for Multi-dimensional 
domain with reflecting boundary. Probability Theory and Related 
Fields,74 (1987)455-477. 

[19] Stroock, D.W. and Varadhan, S.R.S.: Diffusion processes with bound- 
ary conditions, Comm.Pure Appl. Math. 21(1971), 147-225. 

[20] Tanaka,H.: Stochastic differential equations with reflecting condition 
in convex regions. Hiroshima Math. J. 9, (1979)162-177. 

[21] Werner. W.: Some remarks on perturbed Brownian motion, in: 
Seminaire de Probabilities XXIX, Lecture Notes in Mathematics, 
1613(1995)37-43. 

[22] Zambotti, L.: A reflected stochastic heat equation as symmetric dy- 
namics with respect to the 3-d bessel bridge. Journal of Functional 
Analysis, 180(2001)195-209. 

[23] Zhang, X.: Skorohod problem and multivalued stochastic 
evolution equations in Banach spaces. Bull. Sci. math. (2006). 
doi:10.1016/j.bulsci.2006.05.009. 



25 



